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Abstract 

The aim of this article is to derive surface wave models in the presence of surface 
tension and viscosity. Using the Navier-Stokes equations with a free surface, flat 
bottom and surface tension, we derive the viscous 2D Boussinesq system with a 
weak transverse variation. The assumed transverse variation is on a larger scale 
than along the main propagation direction. This Boussinesq system is proved to be 
consistent with the Navier-Stokes equations. This system is only an intermediate 
result that enables us to derive the Kadomtsev-Petviashvili (KP) equation which is 
a 2D generalization of the KdV equation. In addition, we get the ID KdV equation, 
and lastly the Boussinesq equation. All these equations are derived for non-vanishing 
initial conditions. 

Subject Class: 74J15, 35Q30, 76M45, 35Q35 

Keywords: water waves, shallow water, Boussinesq system, viscosity, KdV equation, sur¬ 
face tension, KP 


1 Introduction 

1.1 Motivation 

Understanding the evolution of water waves is a longstanding problem. For instance, 
Stokes already discussed a precise wave at least 130 years ago [18j and his name has been 
given to this wave. Such a wave was still discussed more than one century later |3J ! 

One of the striking phenomena is that there exist surface waves that travel almost 
without modification. Such a behavior is an impetus to investigate the dynamics of these 
waves. 

*The author wants to thank Luc Molinet and Jean-Claude Saut for fruitful discussions. 


1 





The first motivation is to understand how they may either disperse or dissipate and yet 
not vanish. It was discovered that there is a regime in which nonlinearity compensates 
dispersion (with no dissipation). 

Another motivation is more technical. Should a wave keep its shape for a long time 
without vanishing, it could be used in telecommunication without requiring repeaters. 
Such a new technology could drastically reduce the cost of telecommunications. 

Beyond mathematical and industrial considerations, mathematicians and physicists 
model the water waves with a fluid flow in a channel of (often) flat bottom, without 
meniscus (appart in |H]). Since the Euler or Navier-Stokes equations are rather prone 
to faithfully model such a fluid, one could be satisfied with either of these models. Yet, 
the Direct Numerical Simulation is too expensive and asymptotic models are required. In 
such models, one assumes a regime of dimensionless parameters’ smallness and makes ex¬ 
pansions of the equations to derive a simplified model with less fields and less dimensions. 

1.2 Literature 

The techniques are very different depending on whether one chooses Euler or Navier- 
Stokes as a model for the fluid. 

On the one hand, for an inviscid flow (Euler), one assumes very often irrotationality. 
Then one may use a potential function and a Dirichlet to Neumann operator in the 
Zakharov-Craig-Sulem formulation. Numerous articles use this formulation and a review 
of the mathematical proofs has been published HU- Some authors use a velocity-pressure 
formulation, but they still assume the irrotationnality assumption jlj. The derivation of 
Kadomtsev-Petviashvili (KP) assumes a predominant propagation direction along x and 
a weakly transverse propagation along y. It is initially done for an inviscid fluid in |[6j and 
the justification of this approximation is done in [TT] (subsection 7.2). In [TO], D. Lannes 
proves that the sum of one wave propagating to the right and one to the left, both obeying 
a KP equation converges (when e tends to zero) on [0, T 0 /e] to a function consistent with 
the Boussinesq system, meaning that it is a solution, should it be only up to 0(e 2 ). 

Various difficulties must be considered with these surface wave models. Firstly two 
scales of the transverse velocity are choosen in the literature. Either one assumes the 
transverse velocity to be 0(1) or 0(y/e). This problem was raised in Remark 3 of |9j and 
here in Remark [3j It leads to different Boussinesq systems. 

In [15] . the authors prove that a dimensionless water wave system in an infinite strip 
under the influence of gravity and surface tension has a unique solution on [0, T/e]. More 
precisely, if the initial solution is sufficiently regular and small (0(i/e)), there exists 
a solution, on this time interval [0, T/e], that will remain small (0(i/F)). In [TB], the 
same authors prove that on the same time interval, these solutions can be accurately 
approximated by sums of solutions of two decoupled KP equations. It enhances the 
results of [9] by taking surface tension and a variable bottom into consideration. 

The Boussinesq equation is a generalization of the KdV equation for waves moving 
both to the right and to the left. From the Boussinesq system, instead of changing to a 
frame moving to the right as is done for KdV, we remain in the same coordinate system 
and get higher order equation (perturbed wave equation), as is done in [5] (p. 216-219). 
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On the other hand, when one assumes the fluid to be viscous, one may not use the 
potential function. Moreover the number of boundary conditions changes then. Viscosity 
has been considered in water waves since 1895 by Boussinesq pQ. The dynamic of viscous 
water waves on finite depth was investigated more recently by Kakutani and Matsuuchi 
[7]. They derived the viscous KdV equation from vanishing initial conditions but made 
the error of using a Fourier transform in time while the problem is of Cauchy type. 

Liu and Orfila [T3], and the coauthors of P.L.-F. Liu in subsequent articles, investigated 
the viscous Boussinessq system with vanishing initial condition, and validated their KdV 
equation by some experiments. For instance they experimentally proved the reverse flow 
in the boundary layer, which is predicted by theory. 

Although it was not done in Sobolev spaces, the derivation of the viscous Boussinesq 
system (ID and isotropic 2D) and ID KdV equation with a non-vanishing initial condition 
was done in ca without surface tension. It was partially the goal of na not to rely on 
the irrotationnality assumption so as to derive the Boussinesq system and KdV equation 
for a viscous fluid. This could be achieved thanks to the fact that one of the equations 
contains u z = 0(e) which is the zeroth order of the irrotationnality assumption. In [12], 
a more detailed bibliography was given and the reason why all the preceding articles did 
not derive the same KdV equation is explained. 

In this article, we intend to cross the results for a viscous fluid with some results con¬ 
cerning inviscid fluids (KP, Boussinesq equation, surface tension). 

The outline of this article is as follows. In section [21 we present the equations and state 
the viscous Boussinesq system in 2D with surface tension and a predominant propagation 
direction. This result is needed to get the new results on viscous ID KdV equations in 
section [3] the viscous KP equation (2D generalization of KdV with a weak transverse 
propagation) derived in section EB and the viscous Boussinesq equation derived in section 
0 Finally, we justify that our Boussinesq system is consistent with the Navier-Stokes 
equations in Section [6j All these derivations are done for a viscous fluid with surface 
tension and non-vanishing initial conditions. 


2 2D geometry with viscosity and surface tension 

Below, we define the geometry, write the equations, and then make these equations di¬ 
mensionless (subsection 12.IF These equations model a propagation predominantly in the 
x-direction. This will enable us to make a linear analysis with the dispersion relation and 
an asymptotic of the phase velocity in subsection 12.21 A discussion of the relevance of 
viscosity and surface tension will be given in subsection 12.31 We will then use the results 
of P2], where the ID viscous isotropic Boussinesq system and KdV equation are derived, 
so as to state the weakly transverse 2D viscous Boussinesq system (subsection 12.41) . Our 
present derivation includes surface tension, viscosity and a propagation mainly along x 
and weakly along y. Since the already published results are very close, most proofs are 
omitted. 
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2.1 Navier-Stokes equations 

Let u = (u, v, w) be the velocity of a fluid in a 3-D domain Cl — {(x, y, z ) / (x, y) G M 2 , z £ 
(—h,fj(x,y,t))}. So we assume the bottom is flat and the free surface is characterized by 
z = fj(x, y, t ) with fj(x, y, t ) > —h (the bottom does not get dry). Let p be the pressure and 
D[u] the symmetric part of the velocity gradient. The dimensionless domain is illustrated 
in Fig. ED We also denote p the density of the fluid, v the viscosity of the fluid, g the 



Figure 1: The dimensionless domain 


gravitational acceleration, k the unit vertical vector, n the outward unit normal to the 
upper frontier of D, a the surface tension coefficient, and Patm the atmospheric pressure. 
The original system reads: 






div u = 0 


— z/Au + Vp 


-pg k 


(—pi + 2z/D[u] ) .n = -patm,n + 
fjf + ufjx + vfjy — W = 0 
u = 0 


in D, 
in D, 

on 5 = fj(x, t), 

on 5 = fj(x, t), 
on 5 = —h, 


( 1 ) 


where we write the second order tensors and the vectors with bold letters. The differen¬ 
tiated functions are denoted either df/dx, f x , or d x f. The surface tension term includes 
R given by geometrical computations: 

1(1+ ''/jpfa + (1 + '/I)')»., - Zrjxifflxs 

R~ (i + al +i?l ) s/2 

Of course, we need to add an initial condition and conditions at infinity. 

So as to get dimensionless fields and variables, we need to choose a characteristic hor¬ 
izontal length l which is the wavelength along the propagation direction, a characteristic 
vertical length h which is the water’s height, and the amplitude A of the propagating 
perturbation. Moreover, let U, V, W, P be the characteristic horizontal velocity along x, 
horizontal velocity along y, vertical velocity and pressure respectively. Since the zeroth 
order equation is a wave equation of dispersion relation oj = c 0 k x (l + ky/2/kl) + o(ky), 
one may infer that a good scaling for the y direction is if the extra term is of the same 
order of magnitude as e. So k y = O and the interesting scale for y is 0(e -1 / 2 ) more 
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than for the x direction. We may then define: 


c 0 — Vgh, £ — p, P — p-, U — ec 0 , V — ec 0 , 

W = ^c 0 , P = pgA, Re = ^,Bo = ^, 

where Co is the phase velocity. Then, one may make the fields dimensionless and unsealed: 

u = Uu, v = Vv, w = Ww, p = patm — pgz + Pp, fj = Ag, 


and the variables too: 


x — lx, y = e 1/2 ly, z = h(z — 1), t — tl/c 0 . (2) 

We also make the Boussinesq approximation (/3 and e of the same order of magnitude) 
and even (3 = e. Notice that we take the scale used by JD] (V = ec 0 ) and not the one 
taken by [5] (V = £ 3 Pcq). The difference will be highlighted in Remark [31 

With these definitions, the new system with the new fields, variables and outward 
unit normal, still denoted n, writes in the new domain Q t = {(x,y,z), (x,y) G M 2 , z G 
(0,1 + eg(x, y, t))} (see Figured]): 


U t + £UU X + £ 3 / 2 VU y + WU 
V t + £UV X + £ 3 Pw v + WV , 


w t + euw x + e 3 Pvw y + 
u x + Jivy + w z /e = 0 


_£ (. 

Re v 

z Re 
Re 1 


(g — p) n 


l 

Re 


y+^)+P X = 0 

in Q t , 

+ + yf = 0 

in Q t , 

] yy + p) + Vz = 0 

in tt t , 


in Q t , 

)x) ^ Z W X ^ 



(eu y + y/ev x ) 2ev y (v z + y/ew y ) 

u z T w x (v z T y/'e.Wy') 2 

\ V£ / 

1 (1 + £ 4 A xx + e(l + £ 3 V 2 x)Vyy ~ 2 e 4 g x g y g xy 


.11 


Bo (1 + e 3 g 2 x + £ 4 p 2 ) 3/2 

g t + eug x + e 3 / 2 vg y — w/e — 0 


.n 


on z — 1 + eg, 


on z = 1 + eg, 

u = 0 on z = 0. 

(3) 

Simple computations give the outward non-unit normal n = (—eg x , —eg y , 1). 


2.2 Linear theory 

The whole subsection below is a straightforward modification of the case of a viscous fluid 
without surface tension, fully studied in [12] , 
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2.2.1 Dispersion relation 


We are looking for small fields. So we linearize the system (J3]) . The dynamic condition at 
the free boundary can then be stated more explicitely. We get: 


(V’xx T £ ^yy T 'U'zz! £ ') ~P P x ^ 

fs 

v t - + ev yy + v zz /e) + y/Ep y = 0 

W t - (w xx + £Wyy + w zz je ) + p z = 0 

u x + y/evy + w z je = 0 
< (u z + w x )/Re = 0 
(v z + y/ew y )/Re = 0 

(h — P) + j;> e ^/7 = jfo(Vxx + £r lyy) 

r) t - w/e = 0 
u = 0 


in M 2 x [ 0 , 1 ], 

in M 2 x [ 0 , 1 ], 

in M 2 x [ 0 , 1 ], 
in M 2 x [ 0 , 1 ], 
on z — 1 , 
on z — 1 , 

on z — 1, 

on z— 1 , 
on z = 0 . 


(4) 


In order to eliminate the pressure, we differentiate (jTlb with respect to z and fldlh with 
respect to x and compute their difference so as to simplify p xz . Symmetrically, we dif¬ 
ferentiate @ 9 with respect to z and (] 4])3 with respect to y and compute their difference 
(with a yfe coefficient). Pressure is eliminated in two equations on three: 


^ zt W.xt p>g ^xxx) T £ i^yyz ^yyx ) T“ (Uzzz ^xzz)I&) 0) 

Vzt '\/ £r UJyt ppg ((^raz \f £ ^xxy ) T £ {^yyz yf £ '^yyy) "P (Uzzz \f £ '^yzz) j O’ 

So as to eliminate u and v thanks to the incompressibility (HL from the previous system, 
we differentiate the hrst equation with respect to x and the second with respect to y. 
After some simplifications, a good combination of the two gives (thanks to (JUki): 


(d 2 + £&l + e 2 d 2 )(-Rey/ed t + {d 2 z + ed 2 x + e 2 d 2 ))w = 0. (5) 

Let w be of the form A(z) exp i(k x x + k y y — out) with non-negative k x , k y and a (complex) 
pulsation oo. We can define a complex parameter with non-negative real part, similar to 
the one used by [7j: 

p 2 = ek 2 — iuRey/e, ( 6 ) 

where k 2 = k 2 + ek 2 . Thanks to this notation, the solutions of (EJ) are such that 


A{z) = C\ cosh yfek[z — 1 ) + C?, sinh yfek{z — 1 ) 

+ C 3 cosh fi[z — 1) + C 4 sinh p[z — 1). (7) 

The surface tension did not appear yet because it is only in the boundary conditions. 

Up to now we have eliminated u, v and p only in the volumic equations. We still have 
to use the boundary conditions of (JU) to fold the conditions on the remaining held w. 

Easy computations on the boundary conditions of (HJ), similar to the ones done in [12], 
give: 

■uq(O) = 0, 

w( 0 ) = 0 , 

w zz (l) - ew xx (l) - e 2 w yy { 1) = 0, ( 8 ) 

W xx T" lUztti, l) T ppg i^xxzt T" £r ^yyzt) T j^g^,^zzzt (1) j > _ ( k) x T £ dy) W. 
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The solutions 0 satisfies a homogeneous linear system in the constants Ci, C 2 , C3, C4. 
Its matrix is: 


( v^ifesinh (fek) — fek cosh (fek) /i sink /1 —/j cosh /j \ 

cosh (fek) — sinli (fek) cosh/i — sinli // 


V 


2k 2 e 

uz _Li. 

Bo 


0 


feu 2 k + 


2ieujk 3 


/i 2 + ek 2 
-k 2 -^ 


0 

2^y/eiu}k 2 

Re 


(9) 


/ 


Re Bo 

where fc 2 = fc 2 + £& 2 It suffices to compute its determinant to get the dispersion relation: 


4 ek 2 p(ek 2 + p 2 ) + 4 pk 3 e 3 ^ 2 (p sinli (ky/e) sinli /i — cosh (ky/e) cosh /i) 

— ( ek 2 + /f 2 ) 2 (/n cosh [kyfe) cosh — fc-v/isinli {kyfe) sinli //) 

— (fc + fc 3 /Bo)v / £Re 2 (/i sinli {kyfe) cosh /j — kfe cosh ( kfi) sinli /i) = 0. (10) 

This relation is identical to the one of ra appart from the surface tension term and 
the new definition of k for the 3D geometry ( k 2 = k 2 + ek 2 ). Our process of 11011 - 
dimensionnalizing makes a difference between x and z. This is why [7| have k terms, and 
we have kfe terms instead. 


2.2.2 Asymptotic of the phase velocity (very large Re) 

In this subsection, we state the following Proposition concerning the phase velocity: 
Proposition 1. Under the assumptions 

kfe Re c —» + 00 , 
k= 0(1), 
e —y 0, 

Re —> + 00 , 

c = 0(1) (and c bounded away from 0), 

Bo bounded away from 0, 


if there exists a complex phase velocity c = uo/k solution of HI Ch) . then it is such that: 


c = 



k 2 \ tanh {kfe) e m ^ Re 1 ^ 2 {kfe) 1 ^ 


kfe 


2 tanh 3 / 4 (kfe) 


+ o(£" 1/4 Re' 1/2 ). (11) 


Moreover, the decay rate in our finite-depth geometry does not depend on the surface 
tension but only on the viscosity. It is: 


Im(oj) = Im(kc) 


— _ fc5/4gl/8 _+ 0 O- 1/4 R e - 1/2 

2^2 fm tanh 3 / 4 (kfe) { 


( 12 ) 


Notice that our geometry is not infinite. So our formula may be different from Boussi- 
nesq’s and Lamb’s. 

Appart from the surface tension terms and the change of notation of k, the proof is 
not different from the one of Proposition 1 in [12]. 

To what extent is surface tension relevant ? 
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2.3 Relevance of surface tension and viscosity 


As was explained in the previous subsection, surface tension yields only a real term in 
the phase velocity. So it does not give dissipation, but only dispersion. On the contrary, 
viscosity influences both the real and the imaginary part of the phase velocity. 

We choose two different fluids (water and mercury) in order to discuss the relevance of 
surface tension and viscosity. Their physical parameters are listed in Table [Q 



a(Nm 4 ) 

p(kgm 3 ) 

z/(Pa s) 

1/Re 

1/Bo 

water 20°C 
mercury 20 °C 

7.3 10~ 2 

4.4 lO’ 1 

10 3 

1.3 10 4 

10" 3 

1.56 10" 3 

3.2 10 -7 /h 3 / 2 
3.8 10~ s /h 3 / 2 

7.410” 6 /; 2 
3.5 10~ 6 // 2 


Table 1: Physical parameters for two fluids 


Given < 7 , p, u, the Re and Bo number are evaluated in Table [D Depending on h or l. the 
Re or Bo dimensionless parameters are greater or less than 1. The critical values for each 
of them are given in Table [2] It appears that viscosity is not so relevant as surface tension. 



h'crit— Rg 

^crit— Bo 

water 20 °C 
mercury 20 °C 

4.6 10" 5 
1.110" 5 

2.710~ 3 
1.9 10~ 3 


Table 2: Critical values of length (in m) for which 1/Re or 1/Bo is 1. 

Furthermore, so as to compare 1/Bo and e, and under the Boussinesq approximation, one 
may compute the ratio 

J_/ = ° 

Bo £ pgh 2 ’ 

The critical value of h, for which this ratio is 1, is h app = yja/(pg). This other critical 
value is either 2.7 mm for water or 1.9 mm for mercury. In other words if h > h app , then 
the surface tension term is irrelevant and even less important than gravity measured by 
£ in the equation. In most current experiments, h is more than 0.1m and so, surface 
tension should (globally) not be taken into account. As D. Lannes says: “for coastal 
waves of characteristic length L x = 10 m, capillary effects represent only 0.0003% of the 
gravity effects” m section 9.1.2). 

Surface tension is irrelevant for flows with h 0 more than 1 cm. Yet, even if the initial 
flow satisfies this, it may or may not remain true. For instance, in case of flows with 
ripples created by the wind or in case of wave breaking, very short wavelengths appear. 
See section 9.1.2 of m for some references. 

In the case of a flow in a shallow channel (some mm height), and for mercury, a 
depression solitary wave was predicted (as early as the initial article of Korteweg and de 
Vries 0) and observed in [2], In this experiment, A = 0.064 mm, h = 2.12 mm. As a 
consequence, Re=2547 and viscosity should be taken into account. 

In the following, we assume 1/Bo between 0(e) (/ ~ 1mm) and 0(e 2 ) (/ ~ 1cm). 















2.4 Viscous Boussinesq System 

Let us first discuss the regime binding e, Re, Bo. 

The complex phase velocity ffTTl) contains gravitational and viscous terms that we want 
to compare so as to find the regime at which their variations are of the same order of 
magnitude. The first term is the gravitational term (a/( 1 + A; 2 /Bo) tanh (ky/e) / (kyfe)) 
which may be expanded when e tends to zero: 1 — k 2 e( 1 — 3/(eBo))/6 + 0(e 2 ) + 0(Bo^ 2 ). 
The second is purely viscous and can be expanded: — \/2(l + i)(4\/fc) -1 (Re y/e)~ 1 ^ 2 + 
o(Re y/e)~ 1 / 2 ). So if 1/Bo = O(e) or less but l/(eBo) b 1/3, the variations of c on 
the gravitational and on the viscous effects are of the same order of magnitude when 
e(l — 3/(eBo)) and (Re-^/i) -1 / 2 are of the same order. In this regime of very large Re, 
studied hereafter, the dependence of Re on e is such that (to simplify, we assume here 
1/Bo = o(e)): 

Re - £~ b/2 . (13) 

This would be wrong if l/(eBo) ~ 1/3 as is well-known for non-viscous fluids. In that 
case, we would need to go further in the expansion. But since it was assumed above that 
1/Bo is between 0(e 2 ) and 0(e), we are driven to choose the regime (ITHl) and exclude 
l/(eBo) ~ 1/3. 

Our main purpose here is to state an asymptotic system of reduced size from the 
global Navier-Stokes equations in the whole moving domain with surface tension and a 
predominant propagation direction. 

Proposition 2. Let x = (x,y) G M 2 . Let rj(x,t) be the free boundary’s height. Let 
(u b,0 (x, 7 ),u 6,0 (x, 7 )) /or 7 G (0,+ 00 ) (resp. (u u,0 (x, z),v w,0 (x, z)) for z G (0, l + £ 7 ?(x, t))) 
be the initial horizontal velocity in the boundary layer (resp. in the upper part of the 
domain). Let /(x, 7 ) := u b,0 (x, 7) — u u,0 (x, z = 0) and g(x, 7 ) := v b,0 (x, 7 ) —u“’°(x, z = 0). 
If for any xGK 2 , 7 G/ and 7 1 —>- <7 are uniformly continuous in 7 and f x , g y G L* (M + ); 
then the solution of the Navier-Stokes equation with this given initial condition satisfies 
the weakly transverse Boussinesq system: 


u t + 77* - 


+ £UU X + £ 3 / 2 VU. u — 


tt 


[z 2 - 1) 


= 0(£ 2 ) + 0(£/ Bo), 


Vt + V^Vy - + £UV X + £ 3/2 VVy - £r) ytt ^ ij — = 0(£ 2 ) + 0(s/ Bo), 

Vt + U x (x, Z, t ) - | r) xxt (z 2 - |) + e(ur]) x + yj~£v y + £3 (vrj)y - —Ii=(u x + y/£Vy) * -i= 

e -7 ^7^7" = o(e), 


2 £ 


»+oo 


div 


7T 


' 7"=0 


(u b,0 (x, 7") — u u,0 (x, Z = 0)) 
^(yb,°(x, 7 ") _ V u ’°(x, Z = 0)) 


ib" 


1 Y—0 


(14) 

where the convolution, denoted with *, is in time, the parameters £, and R have been 
defined and z G (0,1 + £rj(x, t)). 


With the change of notation k 2 = k 2j r k 2 , the whole proof is the very same as in [12] 
and is omitted. 

Had we proceeded in the same way as [7], we would have distinguished two subdomains: 
the upper part (z > e) where viscosity can be neglected, and the lower part (0 < z < e) 
which is a boundary layer at the bottom and where viscosity must be taken into account. 
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The resolution in each part would have enabled to match the solution on the common 
boundary. Instead, we assume overlapping domains. 

Remark 3. As is stressed at the non-dimensionnalizing step, there are two choices for the 
scale of the transverse horizontal velocity. Either one takes V = £ 3 ^ 2 Cq (see Johnson’s book 

for instance or many others), or one uses V = ec o (see the article of Lannes-Saut fflj 
for instance). Firstly one must notice that since this field is not eliminated (as w is), its 
scale is relevant. Moreover the latter scaling is more general than the former one because 
it does not assume v to be O ( \fe). The main difference is that v Lannes-Saut = V^ v Johnson, 
where one implicitely assumes that the fields VLannes-Saut and vj 0 h nS on are 0(1). To be 
clear, Lannes-Saut’s choice amounts to assume v = 0(1) although it can be written yflidyxf 
! If we plug this v = y/edyif in [lffi , we get the same transverse Boussinesq system as 
1 3J/ (p. 210) because £ 3 ^ 2 VLannes-Saut transforms into e 2 vJohnson and go to 0(e 2 ). So the 
non-dimensionnalizing process of Lannes-Saut is weaker than the one done by many other 
authors. 

Lannes and Saut m give one more argument by exhibiting an (odd) solution to the 
linearized system (in their Remark 3) such that v = y/Edyif. It is of size 0(1) and not 

0(y/£): 


u = 0, v = f\y - y/et) - f\y - y/et),^ = f\y - y/et) + f\y - y/et). 


We want to notice here that the potential representation of the velocity might have induced 
the idea that if the potential is physical, then the deduced field (the velocity) is such that 
(u,v) = (d x (j), y/edyfi) and so, that v = 0(y/e). Had we started from the velocity pressure 
formulation, this assumption would not be natural. Yet, such an attempt of explanation 
does not match reality. Johnson /3/ uses the velocity pression representation as we do and 
assumes v = 0(y/e), while Lannes and Saut use the potential representation and do 
not make this assumption ! 


Remark 4. It was already noticed in fWj that the double integral term in (If) is new and 
surprising because of its dependence on the initial condition. If we assume an initial flow 
of Euler type (so that u b, ° — u u, ° = 0), the double integral vanishes. But is it physical ? 
In other words, does an initial (inviscid) flow in the boundary layer (where Navier-Stokes 
applies) establish (as a Navier-Stokes flow) fast or not ? Can we assume an initial Euler 
flow and a Navier-Stokes evolution PDE without loss of generality 1 

We claim that the answer is negative for at least two reasons. On the one hand, the 
characteristic time for the viscous effects to appear is roughly T^se = ph^/v or T^se = 
pi 2 /u. Then, its ratio with the characteristic time of the inviscid gravity flow (1/cq) is 
either Re yfe = e~ 2 or Re/y^e = respectively. Whatever the chosen characteristic 
time, this ratio is large and the flow in the boundary layer does not establish fast enough. 

On the other hand, as is argued in m. the value of this integral for a typical exponential 
flow in the boundary layer can be computed. It tends to zero only like 1 /y/t as is classical 
for Navier-Stokes flows. As a consequence, one may assume that it goes from zero (Euler) 
to exponential (boundary layer) within the time (1 /y/t) which is large with respect to the 
characteristic time of the flow. 
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3 Viscous KdV with surface tension 


Starting from fj!4f) . one may assume the fields do not depend on y , and v — 0. So the 
flow is purely one dimensional (one direction). One may see that the zeroth order of 
this Boussinesq system is the wave equation. Then one knows that there are two waves 
propagating in each direction. If we look only for the waves that propagate to the right, 
one may make a change of variables suggested by the zeroth order equation: 


£ = x — t, t — et. 


Every term can easily be converted in these new coordinates appart from the convolution 
(one derivative and a half integration), the new surface tension term and the double 
integral. The first was treated in a very clean way in [12]. The second is new, but very 
simple. The third one was treated in H 2 , but its treatment is improved below. After a 
change of variable, this term writes: 


+ ~j= [ ( u ®°(£ + ~i l") ~ m x’°(£ + ~,z — 0)) x f 

\/7T ./y/_o V £ £ / ./„/ 


'Re it 
4 T ' 


,~T 


V/dV' 


(15) 


Y=o 


The article H2| argued that if the initial horizontal velocity is localized and r not too 
small, then u b,0 (f + r/e, 7 ") — u u,0 (f + r/e, z = 0 ) will be negligible in comparison with 
e. More precisely, if £ (-»• w 6,0 (£, 7 ") — u u,0 (f, z = 0) tends to zero when £ tends to +00 at 
least as £ _1 . Then the term (1131) is at least 0(e 2 ). 

One may add one more argument if r is still not too small. The upper bound of the 
inner integral contains a yfe. So mainly large 7 " (larger than 1 j\fe) will be relevant. But 
for those 7 ", the w fc,0 (£ + r/e, 7 ") — u u,0 (£ + r/e, z = 0 ) is small uniformly in £ + r/e 
because of the matching condition on u° between the boundary layer and the upper part. 
So (1131) is small for two different reasons and we can state the viscous KdV equation with 
surface tension in the following Proposition. 


Proposition 5. If the initial flow is localized, the KdV change of variables applied to the 
ID version of the system [Iff) leads to 


2r] T + + 




1 r /e %«+M <=0(1)i 


3 eBo J yj 7 tR j£i = 0 vV 

for not too small times r, where we set Re= Re~ 5 l 2 . 


(16) 


We do not replace the upper bound of the convolution {r/e) by +00 for two reasons. 
On the one hand we do not know how fast this integral converges when r/e tends to + 00 . 
On the other hand the r term reminds us that this integral on £' is a mixture of time and 
space. 


4 Viscous KP equation with surface tension 

Below, we derive the viscous KP equation, then discuss the zero-mass (in £) constraint. 
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4.1 Derivation 

We recall that we set Re= Re ~ 5//2 , and 1/Bo = 0(e) but l/(eBo) ^ 1/3 + o(l). 
With these assumptions, the change of variables 


f = x - t, y = y, r = et , 


and the assumption that the initial flow is localized, the system (fDl) writes (since e/Bo = 
0 ( 5 2 )): 

-uz(€, y, Z, r) + + e(u T + uu^ - r)^(.z 2 - l)/2 - ^ + y/evu y ) = 0(e 2 ), 

-V£ + \[er)y + e(v T + uv£ - yjer] y ^(z 2 - l )/2 - 7 ^- + \[evv y ) = 0 (e 2 ), 

y eHo 

+ £(Vt + - V 3 ) + ( w l)n + Ve(vri)y) + Vev y ^ ^ 

g [ T/£ uflf + f,',y,z,T) + Veu y (£ + £',2/,z,T) , = Q 2 

VnR j £=0 V^ 7 

In the above system, we did not write the double integral term from ffT5|) . Indeed, it was 
justified in the previous section that this term can be neglected. Moreover, one might 
prove that the transverse velocity is such that v(z) = v(z') + 0 (e) in a similar way to 
Lemma 11 of [T2| (roughly, one differentiate (1 14 H ? with respect to z, then integrate with 
respect to time t). So as to go further, since, at the first order u% — r/y one may justify 
that u T = rj T + o(l) as in the inviscid case, and u — ij + o(l). One may then add the hrst 
and third equation: 


2er ] T + 3er/r/ ? + 77 ^ + e 3/2 vu y + e 3/2 (vrj) y + y/ev y 

_ [ T/e Vt(Z + r ) + y/evyit, + V, z » r ) dr / = q, 2s (18) 

v^rJ ?=0 ve ? 

In order to eliminate v, one must differentiate with respect to £ to use the corresponding 
(TTTlB . Our choice of non-dimensionalization forces us to manage extra terms, but we get 
the classical KP equation with a viscous term: 


2 r] T + 37777 ^ + 


eBo 


r te5 +1 lyy- 


1 r /£ + 

yfxR J £'=0 V? 


= -y/e (2 vrj iy + + n Ty ) + 0(e) + 0( 1/Bo), 


where n satishes (ll7P o. This completes the proof of the following Proposition. 


Proposition 6. If the initial flow is localized, in the Boussinesq approximation and if 
1/Bo = 0(e) but 1/Bo 7 ^ e/3 + o(e), the surface waves of a viscous fluid propagating 
predominantly along x and weakly along y satisfy the viscous KP equation: 


(^Vt + 3777fc + 




+ Vyy 






O(Ve), (19) 
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or the equivalent mixed system 


2i) T + 3 r)r)£ + 




~ 




vd£ + £',y, T ) 


dF + Vy 


-V£ + r)y 


0(y/e), 

0(\fe). 


( 20 ) 


One must keep in mind that our choice of scale triggers the yfe (even inside the con¬ 
volution) in (1171) . This choice is already discussed in Remark [3J Had we made the same 
choice of scale for v as Johnson [5] and many others, we would have replaced v by y/evj 
in (1171) . The equivalent (1T81) could be simplified by £ and it would write : 


2 Tjr + 37777$ + 


, 1 f T/e vt{£ + ?,y,T) N 

+ v Jv - -J= J ^^-df = O(e), 


3 eBo 

and vj would satisfy 1132 modified: 

-Vj£, + Vy + £ (vjt + uvj£ - Vy^(z 2 - l)/2 - -^0 = 0(s 3/2 ). 


Oddly, this would lead to the same KP equation as (H9l) but up to the order 0(e) and not 
0(^/e). During the proof, we exhibit the 0(yfe) terms which would be 0(e) if v = \fevj 
although they are not kept in the Proposition. So there is no contradiction and the 
accuracy of KP is tied to the property that v = 0(y/e) or not. 


4.2 The zero-mass constraint 

As is very well discussed in im the usual KP equation is often written with an operator 
d^d 2 . Yet such an operator assumes the solution is differentiated (in £) from a function 
that tends to 0 when £ —>■ ± 00 . It is proved in da that although it is not obvious, this 
assumption is right for inviscid fluids. Their proof is done in two steps. 

In the first step, the linear KP equation is solved with inverse Fourier transforms. For 
more general equations, but only dispersive, like u t — Lu $ + d^ l d yy u = 0 (the symbol of 
L is e | X | Q for a > 1/2), the fundamental solution writes 

G(t,£,y) = ^X,y>)^,y) 

where e = ±1 depends on the KPI or KPII equation (or on the sign of 1/3 — l/(e:Bo)). 
For the usual KPI and KPII, a = 2. 

Thanks to the Lebesgue’s Dominated convergence theorem and inventive changes of vari¬ 
ables, it is proved that the fundamental solution is regular, and is differentiated from a 
more regular function A (which is C(M 2 ) f) L°°(M 2 ) P|C|(M 2 )) such that A * u 0 —» 0 when 
| £ |-» Too. 

In a second step, the nonlinearity is treated with the Duhamel formulation. The 
solution is then the convolution of the group of the linear equation with the nonlinear 
term u and all what is already proved for the linear equation extends to the nonlinear 
one. 


e it(ex\x\<*- y 'yx) 
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One may wonder whether such an interesting result can extend to the viscous case. 
Indeed, the proof is even simpler and the details are left to the interested reader. The 
Fourier transform of the dissipative term writes 


^Xr' F ^ x 



vt{€ + ?,y,T) 

VF 




d£'d£ 


By a simple change of variable, one may see that the last integral is indeed only a function 
of the sign of X (denoted sgn(A)) times 1 /y/1 X |. This function of sgn(A) can be 
computed : 

r+°° P i sgn(x)£ nr 

L 2 {i+isga{x)) ’ 

and the Fourier corresponding term reads: 

-j====(l +isgn(A)). 

As a consequence, one could study the fundamental solution 


G NS (t,^y) 




' e it(sX\X\°‘-y' 2 /X+(l+isgn(X))X(2R\X\)- 1 -/ 2 ) 


The only real part (non-dispersive) inside the exponential is it x fsgn(A)A(2i? | X \ 
)^ 1//2 — —t \ X I 1 / 2 (2 R)~ l t 2 which sign is compatible with the dissipation and ensures 
convergence of all the integrals. So one do not even need to use the Lebesgue dominated 
convergence theorem and the proof is simpler than in HU. 


5 Viscous Boussinesq equation with surface tension 

The Boussinesq equation is a second order equation that takes into account waves going 
both to the right and to the left. Since the derivation is straightforward, we only sketch 
it, following Johnson |5j (p. 216-219) and state the Proposition [3 Johnson [5j uses a 
different scaling, but only along the transverse direction which is not used here. 

Starting from the Boussinesq system (fT4l) written only in ID (with x,z,t and not 
x,y,z,t ), one may differentiate (HDD with respect to t and (fTDh with respect to x. The 
difference of these two equations writes after some easy computations: 


Vtt VxX 


£ I U 2 + 


£ TJxxxx £ 1 

T^Vxxtt H ^ I ,—— Vtt * ~~F 

3 Bo yfnR y/t 


+ 


£ 1 
=p.v. 


VtiR Vi 


r»+00 


u x’°( x , z = 0) — / u b f(x,V)e-^ V 


'7'=0 


= 0(e 2 ), (21) 
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where the convolution is in time and p.v. denotes the principal value as defined in the 
theory of distributions. The p.v. stems from an integration by parts in which one has a 
boundary integral 


£ 

\J 7T Rt 


/°(u7') e ~ 


i?7 z 
4 1 


“1 +oo 


7'=0 


The u u, ° and u b, ° are the initial values of u u and u b . They must be provided. Ap- 
part from the two viscous terms and the term generated by surface tension, this equa¬ 
tion is identical to (3.41) of [5]. Since we concentrate on x = 0(1), we may write 
u = n 3 UX J-oc 5 rj t (x', t)dx' + 0(e). The above equation is then the Eulerian form of 
the Boussinesq equation rewritten in (1^1) . 


Let us now come back to the Lagrangian coordinates and change of fields. In [5], the 
author proposes X = x+e ff rj(x', t)dx', but we consider the following change of variable 
easier to justify: 

X = x-ef*u(x,t')dt' 

H(X,t ) = H(x — e f* u(x, t')dt', t) — ri(x, t) — er] 2 (x, t). 

The function H is assumed to be dehned on X in the initial domain. Using these defi¬ 
nitions, it is straightforward to rewrite (12TT) under the form of a Lagrangian Boussinesq 
equation with surface tension and viscosity taken into account. Only the term containing 
the initial horizontal velocities needs some insight, ft is the last one in (H?T1) and both 
u u ’°(x,z = 0) and u b ^(x,fl) can be considered as depending on x or on X since an £ 
appears in front of the whole term and the difference between x and A" is 0(e). The 
following Proposition states the two forms of the Boussinesq equation in ID. 

Proposition 7. Let a flow of a viscous fluid with surface tension in a 2D channel with a 
flat bottom like the one depicted in FigureUi Let u u ’°(x, z = 0) denote the initial horizontal 
velocity in the upper part (respectively u b,0 (x, 7 ) in the boundary layer). We assume the 
Boussinesq approximation and 1/Bo = 0(e). Under these assumptions, the surface waves 
obey the viscous Boussinesq equation in the Eulerian form: 


hit Qxx ^ 



Vx 


+ 


V£r , u * Vi 


+ 


WR^Vi 



(Vt)t =0 - 




0(£ 2 ), (23) 


where the convolution is in time and p.v. denotes the principal value as defined in the 
theory of distributions. The Lagrangian form (with f22 [)/ writes: 


3 


Htt(X,t) — H \x — -e(H 2 )xx — ( 7 — — ) Hxxxx + 


1 

Bo 


-=//« * -1= 
\ZttR \/t 


e 1 

+ ^=p.v.— 
V nR y/t 


-(H t (X,t)) t=0 - 


^+00 


u 


b f(x,fl) e-ih' 2 d 7 ' 


*7 ,= 0 


= 0(£ 2 ), (24) 


where the initial velocities defined on x can be evaluated either on x or on X. 

Notice that the formulae (I23II24D also have a formulation in terms of u u ’°(x, z = 0) 
u b,0 (x, 7 ) which is physically meaningful. 
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6 Is the Boussinesq system consistent with Navier-Stokes 

? 


The weakly transverse Boussinesq system is stated in Proposition [2J We intend to justify 
here that it is consistent with the initial Navier-Stokes equations (J3J) . 

What is consistency ? In na, the author justifies, but does not prove in any functionnal 
space, that if there is a solution to Navier-Stokes equations, with a free boundary and a 
flat bottom, then, under the Bousinesq approximation, the solution satisfies the viscous 
isotropic Boussinesq system. Such a necessary result deserves to be completed by a 
sufficient one which is consistency. This is the goal of the present section. 

Notice that the consistency as defined by Lannes in [IT] (Definition 5.1) is much more 
rigorous. While we use only expansions, Lannes asked whether a solution to the asymp¬ 
totic model (assuming it exists), in the asymptotic regime, satisfied the initial complete 
model on a convenient time-interval and in convenient norms, up to a given power of the 
small parameter. Lannes proved in [10] that KP solution (with his scaling) is consistent 
with the Boussinesq system. 

We prove below the following Proposition in the ID case. 


Proposition 8 . Let u u (x,z,t) be the horizontal velocity in the upper part of a flow in a 
geometry as in Figure [H and 77 (x, t) be the height of the free boundary above the fluid, for 
x = (x,y) € M 2 ,z e (e, 1 + eg(x,t)) and nonnegative time t. Let (u u,0 (x, z), v b,0 (x, z)) 
denote the initial horizontal velocity in the upper part (respectively (u b,0 (x, 7 ), v b,0 (x, 7 )) 
in the boundary layer with 7 = z/e). 

If (u u , v u , 77 ) satisfies the non-isotropic Boussinesq system a then there exist fields in 
the upper part (u u ,v u ,w u ,7 7 ) (resp. (u b ,v b ,w b ,7])) satisfying the Navier-Stokes equations 
0 in the upper part (resp. the boundary layer) of the domain and the interface continuity 
at least as 0(e). 


Proof. We start from the final Boussinesq model (TT4l) written in ID: 

Vxxx (' 


< + Vx 


Bo 


+ £U u U u x - £71 xtt - 


^ =0(e 2 ) + 0(e/ Bo), 


Vt + u x (x, z , t) - ^g xxt (z 2 ~ ~) + e(u u g) x - 


2e 


"+00 


7 r 


'7"=0 


{u b f(x,fl')-u u f(x,z = 0 )) 


V^R Ux * Vt 


IV' 


e ^dq'dq" = o(e), 


1 y=o 


(25) 

where u u is the horizontal velocity in the upper part (z G (e, 1 + eg(x,t))), g(x,t) is the 
free boundary’s height, u u,0 (x, z) (resp. u b,0 (x,^)) is the initial horizontal velocity in the 
upper part (resp. in the boundary layer). In the boundary layer, we set z = £ 7 . 

First, one must justify that g tt = g xx + 0(e), which is easy from the zeroth order of 
(1251) i. Then, the equation (1251) 1 enables to prove the following Lemma. 


Lemma 9. A localized solution of / 1251) i is such that 

fo 1 U = u(x, z, t) - eg xt (z 2 - l/3)/2 + 0(e 2 ), 
u(x, 0, t) = u(x, z, t) — eg xt z 2 /2 + 0(e 2 ), 
u(x, 1, t) = u(x , z, t) + eg xt (l — z 2 )/2 + 0(e 2 ). 
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The way to prove this Lemma is identical to the proof of Lemma 11 of m One 
differentiates (l25]b with respect to z, integrate with respect to time t, prove that the 
constants of integration vanish if the wave is localized, and use that rj tt = r] xx + 0(e) to 
get u z = er] xt z + 0(e 2 ). Completing the proof is then easy. 

Thanks to the previous Lemma, if we define 


p u (x, z, t ) = (rj - r] xx /Bo) - er} tt (z - 1 ) + £ u xt (x, z", t)dz" d z' 


n.U 

)l+eri 


w a (x, z, t) = e(r] t + er] x u u (z = 1 + erj)) — e J* 

then ( u u ,w u ,p u ) satisfies © 1 , 3 , 5 ,6 up to 0(e 2 ) + 0(ej Bo), and (OIL exactly in the upper 
part. 

We must now check ((3|) in the boundary layer, where we define z = £7 and 


u 


’(x,j,t) = u u (x,z = 0 ,t) + /o + °° fo ( x » i) - ~ d ^ 

-u u (x,0, .) * - & / 0 + °° f 0 (x, i) e -dV, 


where fo(x, 7 ) = u b,0 (x, 7 ) — u u,0 (x,z = 0). It is proved in Lemma 6 of | 12 | that this 
function satishes 


(u b - u u (z = 0))t - (u b - u u (z = 0)) 77 /i? = 0, 

=► u\ J ri) x -u\ 1 lR-r] xxx l Bo = 0(e) + 0(e/ Bo). 


This u h also satishes the initial conditions and the boundary conditions: 


u b (x, 7 , t = 0 ) = u b,0 (x, 7 ), 
u b (x , 7 = 0 , t) = 0 , 

u b (x,^f + 00 , t) = u u (x,z = 0 , t) (continuity condition). 


The function u b is dehned as u h — u u (z = 0) being the solution to the heat equation. 
Because it is only the zeroth order of the conservation of momentum in the boundary 
layer @ 1 , and also because the limit between the boundary layer and the upper part is 
not made sufficiently precise through the lift “function” u u (x, z = 0, t), one may not prove 
more. 

If we define: 

P b (x , 7 , t) = rj(x, t ) - r) xx /Bo + er] xx /2 
w b (x, 7 , t) = -e 2 f Q 7 u b (x, y, t) dy, 

then P)i in the boundary layer reduces to: 

(u b - U u )t - (u b - U u )jj/R- eu u u x + eu b u b x - eify fj u b + £rj xxx /2 + 0(£ 2 ) + 0(ej Bo), 
= -eu u u x + eu b u b x - eu b fj u b + £p xxx l2 + 0(e 2 ) + 0(e/Bo) = O(e), 

where u u is evaluated at z — 0. So (u b ,w b ,p b ) satishes (EJ) 1 only up to 0(e). So as to go 
further, one should use correctors, solve a nonlinear heat equation, make a more precise 
study at the interface between the two subdomains and use functionnal spaces. Simple 
computations prove that (u b ,w b ,p b ) satishes flS]) 3 , 4,7 up to 0(e 2 ) + 0(e) Bo). 

The overlapping of the boundary layer and the upper part makes that consistency 
cannot be justihed more precisely than 0(e). Yet, since the interface condition (for 
instance at z = \fe) between these two domains is satished, the proof is complete, up to 
order 0(e). 

□ 
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7 Conclusion 


We modeled the water waves flow by the Navier-Stokes equations in a 3D geometry with 
a flat bottom and free surface. We assumed the order of magnitude of the velocity in the 
transverse direction and deduced the scaling in this coordinate and discussed its influence 
on the results (Remark [3]). The linear theory gave a new dispersion relation and a new 
phase velocity (Proposition [T]) . Surface tension was discussed and seemed often more 
relevant than viscosity, although these two parameters act very differently. We stated the 
associated Boussinesq system. This intermediate system enabled us to derive the KdV 
equation (Proposition 0) and the KP equation or system (Proposition [H]) • We justified why 
the zero-mass constraint does not raise any trouble for viscous KP, using m. Using the 
previous computations, we derived the Boussinesq equation (Proposition [7]). Reciprocally, 
we proved that, should they exist, the solutions to the non-isotropic Boussinesq system are 
consistent with the Navier-Stokes equations (Proposition[S]). These systems and equations 
are derived for a viscous fluid and take the surface tension into account. 
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